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Abstract 

Consistent perturbation theory for thermodynamical quantities in strongly 
type II superconductors in magnetic field at low temperatures is developed. 
It is complementary to the existing expansion valid at high temperatures. 
Magnetization and specific heat are calculated to two loop order and compare 
well with existing Monte Carlo simulations, other theories and experiments. 
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Typeset using REVT^ 



I. INTRODUCTION 



Thermal fluctuations play much larger role in high Tc superconductors than in the low 
temperature ones because the Ginzburg parameter Gi characterizing fluctuations is much 
larger. In the presence of magnetic fleld the importance of fluctuations in high super- 
conductors is further enhanced. First of all, strong magnetic fleld effectively conflnes long 
wavelength fluctuations in direction perpendicular to the fleld so that the dimensionality 
of the fluctuations are reduced by two Moreover highly anisotropic BSCCO type lay- 
ered superconductors are basically two dimensional and one expects further increase in the 
importance of fluctuations. Under these circumstances fluctuations influence various phys- 
ical properties and even lead to new observable qualitative phenomena like vortex lattice 
melting into vortex liquid far below the mean fleld phase transition line. It is quite straight- 
forward to systematically account for the fluctuations effect on magnetization, speciflc heat 
or conductivity perturbatively above the mean fleld transition line using Ginzburg - Landau 
description . However it proved to be extremely difficult to develop a quantitative theory 
in the most interesting region below this line, even neglecting fluctuation of magnetic fleld 
and within the lowest Landau level (LLL) approximation. 

To approach the region below the mean fleld transition line T < Tmf{H) Thouless 
proposed a perturbative approach around homogeneous (liquid) state in which all the 
"bubble" diagrams are resummed. The series provides accurate results at high temperatures, 
but become inapplicable for LLL dimensionless temperature ~ {T — Tmf{H)) / (T HY^"^ in 
2D smaller than 2 (see the lines H4-6 on Fig. 3 for the case oiD = 2 which represent successive 
approximants (the corresponding three dimensional plots appear in |^). Generally, attempts 
to extend the theory to lower temperature by Fade extrapolation were not successful 

Alternatively, more direct approach to low temperature fluctuation physics would be to 
start from the Abrikosov solution at zero temperature and then take into account perturba- 
tively deviations from this inhomogeneous solution. Experimentally it is reasonable since, 
for example, speciflc heat at low temperatures is a smooth function and the fluctuation con- 
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tribution is quite small. This contrasts sharply with theoretical expectations. Long time ago 
Eilenberger calculated spectrum of harmonic excitations of the triangular vortex lattice . 
Subsequently Maki and Takayama noted that the gapless mode is softer than the usual 
Goldstone mode expected as a result of spontaneous breaking of translational invariance. 
The propagator for the phase excitations behaves as 1/(A;^)^ in 2D. This unexpected addi- 
tional softness not only enhances the contribution of fluctuations at leading order but also 
leads to disastrous infrared divergencies apparently at higher orders. As a result, whether 
the perturbation theory around the vortex state is reliable becomes doubtful. For example 
the contributions to free energy depicted on Fig. 2a and 2d are respectively and log^(L) 
divergent {L being an IR cutoff) in 3D and the divergencies get worse at higher orders . In 
2D the situation is aggravated: the diagrams diverge as and L^. Also qualitatively the 
lower critical dimensionality for melting of the Abrikosov lattice is = 3 and consequently 
infinite range vortex lattice in clean materials exists only at T = 0. One therefore tends to 
think that nonperturbative effects are so important that such a perturbation theory should 
be abandoned [Q. However a closer look at the diagrams like Fig. 2a (see details below) 
reveals that in fact one encounters actually only milder divergencies in both 2D and 3D. 
This makes the divergencies similar to the so called spurious divergencies in the theory of 
critical phenomena with broken continuous symmetry. In such case one can prove § that 
they exactly cancel at each order provided we are calculating a symmetric quantity. 

In this paper we show that all the IR divergencies in free energy or other quantities in- 
variant under translations cancel to the two loop order. The 3D case was briefly reported in 

here we present in detail the more complicated 2D case. We calculate magnetization and 
specific heat to this order and compare the results with existing high temperature expan- 



sion, Monte Carlo (MC) simulation of the same system and experiments. Qualitatively 
physics of fluctuating 3D GL model in magnetic field turns out to be similar to that of 2D 
spin systems (or scalar fields) possessing a continuous symmetry. In particular, although 
within perturbation theory the ordered phase (solid) exists only at T = in thermodynamic 
limit, at low temperatures solid with powerwise decay of the translational order or liquid 



(exponential decay) differs very little from solid in most aspects. Therefore, one can effec- 
tively use properly modified perturbation theory to study quantitatively various properties 
of the vortex liquid phase. Similarly physics of the 2D GL model is analogous to that of the 
ID scalar theory with, say 0(2) symmetry. This is equivalent to the anharmonic oscillator 
in quantum mechanics or the XY spin chain in statistical physics. 

The paper is organized as follows. In section II the model is introduced and the Feynman 
rules for the loop (low temperature) expansion within the LLL approximation are set up. The 
free energy to two loops is calculated using diagrammatic expansion in section III. It is shown 
how all the infrared divergencies cancel exactly at least up to the two loop order. In section 
IV we compare the present expansion with the high temperature series, theory of Tesanovic 
et al IjTlJ, available MC simulations and experiments. In section V we discuss qualitatively 
nonperturbative effects using the analogy with magnetic systems mentioned above. We 
exemplify this by calculating perturbatively the ground state energy of the 0(2) invariant 
ID chain (equivalent to the quantum mechanical anharmonic oscillator) which exhibits 
similar IR divergencies. We argue that although the infinite range translational order is not 
present in this system, locally the system looks like a lattice and the perturbative results 
are valid up to exponential corrections. Finally we briefly summarize the results and discuss 
the melting transition observed in experiments and some MC simulations, contributions of 
higher LL, and fluctuations of magnetic field in section VII. 



II. MODEL, MEAN FIELD SOLUTION AND THE PERTURBATION THEORY 



A. 2D Ginzburg - Landau model 



Our starting point is the GL free energy: 



F = L, I — 
2m 



(1) 



Here A = {—By, 0) describes a nonfiuctuating constant magnetic field in Landau gauge. In 
strongly type II superconductors (k ~ 100), far from Hd (this is the range of interest in 



this paper) magnetic field is homogeneous to a high degree due to superposition from many 
vortices. For simphcity we assume a = —a'Tc{l — t), t = T /T^ although this dependence 
can be easily modified to better describe the experimental coherence length. 

Throughout most of the paper the following units will be used. Unit of length is ^ = 
{2mabCi'Tc) and unit of magnetic field is Hc2, so that the dimensionless magnetic field 
is b = B/Hc2- The dimensionless free energy in these units is (the order parameter field is 
rescaled as ijj'^ ^^■?/'^): 



I UJ 



The dimensionless coefficient is 



(2) 



UJ 



V2Gi7rH, (3) 



where the Ginzburg number is defined by Gi = l( 327re ^^€Jc7 y g^j-^^ ^ ^ nic/mab is an 
anisotropy parameter. The coefficient u determines the strength of fluctuations, but is 
irrelevant as far as mean field solutions are concerned. 



B. Mean field solution near H, 



c2 



Define operator Ti. = M—D"^ — b), the free energy becomes 



F _ 1 
T ~ ^ 



d X 



(4) 



Here, 



1-t-b 



ah = 



(5) 



is the second expansion parameter. If at is sufficiently small, GL equations can be solved 
perturbatively (see ll^Jll] for details): 




^(x) + 0(af) 



(6) 



where Pa = 1-16 and 
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y?(x) 



27r6 



exp < t 



7r/(/-l) 27r, /- 

2 OA 



-(v^y 0^ 



OA 



(7) 



The lattice spacing is given by with oa = ^i^iid ^ the magnetic length in our units. 
It is normalized to the unit average < >= 1 

C. Fluctuations and Feynman rules 

Within the LLL approximation, which will be used here, the above solution becomes ex- 
act so higher orders in ah do not appear and ip can be expanded in a basis of quasimomentum 
k eigenfunctions 



27rfe 



around the mean field solution 



2 OA V6 



2 V6 OA J 



V5k(a;)ck (Ok + iAk) 



(8) 



The shift f of the fiuctuating field will be generally different from its mean field value. 
The integration is over Brillouin zone of the hexagonal lattice which has an area 27i. Ok 
and Ak are "real" fields satisfying = 0_k, ^k ~ ^-k- They are somewhat analogous 
to the acoustic and optical phonons in usual solids with some peculiarities due to strong 
magnetic field |^. For example the A mode corresponds to shear of the two dimensional 
lattice. After we substitute eq.(||) into the free energy, the quadratic terms in fields define 
propagators, while the cubic and quartic terms give rise to interactions. The phase factors 
Ck = \J^^^\ with 7k = j^Lpl.{x)Lp*_]^{x)Lp{x)'^{x) are introduced in order to diagonalize the 
resulting quadratic part of the free energy 



quad „ I /o \2 

^ .Ib.z. [2,-n) 



[Po\\^)OlO^ + P^\\^)AlA^]. 



(9) 



Here, Po,A(k) are the propagators entering the Feynman diagrams in Fig. la-b and are 
given by 
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Po,A(k) = Mo,i(k); Mof^(k) = -a, + T;^(2/3k±|7k|), (10) 
For convenience, let us introduce the function A(ki, k2): 



A(ki,k2) = W — cxp<^ } (11) 



2 ' t 2 

V(-)'- exp L— [lk\ + mkl] - \ [(fcf + —If + {kl + — m)^! ) . (12) 

l,m ^ '- 

Now, various functions encountered so far as well as all the three and four leg vertices can 
be expressed in terms of it: for example, 7k = A(— k, k), /3k = A(0, k). The vertices are 
depicted on Fig.lc-i with following expressions. Three leg vertices: 

(c) : AkiAk2A_ki-k2 = -2T;Im [A(k2, ki + k2)Ck^Ck2Cki+k2] , (13) 

{d) : ^ki^kaC-ki-kz = 2vRe [-A(ki, -k2)ckiCk2Ck^+k2 + 2A(-ki - k2, -k2)Ck^Ck2Cki+k2] , 

(14) 

(e) : OkiOk2^-ki-k2 = -2^;Im [-A(ki, 2ki + k2)ck^Ck2Cki+k2 + 2A(ki, -k2)ckiCk2Ckj+k2] > 

(15) 

(/) : OkiOk20_ki-k2 = 2^; Re [A(k2, ki + k2)ck^Ck2Cki+k2] ■ (16) 
Four leg vertices 

{g) : AkiAk2Ak3A_ki-k2-k3 = ^ Re[A(ki + kg, ks + k3)ck,Ck2C^3Cki+k2+k3], (17) 

{h) : AkiAkaOkgO-ki-ka-ks = Re [-A(ki + kg, k2 + k3)ck,Ck2C^3Cki+k2+k3 (18) 

+2A(ki + kg, ki + k2)ckiCk2C*^Cki+k2+k3], (19) 



{i) : OkiC'k2C'k3C'_ki-k2-k3 = \ Re[A(ki + kg, k2 + k3)ck^Ck2C Ck^+ka+ks]- (20) 
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III. CANCELLATION OF INFRARED DIVERGENCIES IN LOOP EXPANSION 



A. One loop free energy, "field shift" and destruction of the infinite range 

translation order by fluctuations 

If the fluctuations were absent the expectation value Vq = ^ would minimize Go = 
—tthv'^ + ^Pav'^- The one loop contribution to the free energy is 

= ^(^/{log[^o(k)] +log[Ml(k)]} (21) 

To this order the free energy which is a symmetric quantity is convergent. However the 
expectation value of the field which is not a symmetric quantity is divergent. Minimizing 
Go + Gi with respect to v would lead to the following correction to the vacuum expectation 
value: 

2_i 1 f / [2/jic+i7ki] , m-bk\] \ ...s 

'^-2{2nyi\ MS(k) + Ml(k) /• ^''^ 

Due to additional softness of the A mode the above integral diverges in the infrared 
region. This means that the fluctuations destroy the inhomogeneous ground state, namely 
the state with lowest energy is a homogeneous liquid in accord with the fact that the lower 
critical dimension for the present model is D = 3 . It however does not necessarily means 
that perturbation theory starting from an ordered ground state is inapplicable. The way 
to proceed in such situations have been found while considering simpler models like ID ip^ 



model F = i(Vv3a)^ + ViVa), a = 1, 2 discussed in detail in section V (see also [P 



B. Two loop contributions to the free energy 

To the two loop order one gets several classes of diagrams (see Fig. 2): the setting- 
sun {AAA,AAO,AOO,000), double bubble {AA,AO,00), and the "correction term" 
{AA, AA,00), which arises due to correction in the value of v from eq.(p^) . All of them 
can be expressed explicitly in terms of the function A(ki, k2) and c/k = 2c|. 
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1. Setting sun diagrams. The setting-sun diagrams are shown in Fig.2a-d. The AAA 
diagram is naively the most divergent one among them. 

^ = i i Ml(p)Mltq)Ml(p + q) ' ^^^^ 

IssAAa{p, q) = |A(p, -q)|^ - X^{p, -q)(ipciqci*+q + 

2A(p, -q)A(p, -p - q)ciqcip+q - 2A(p, -q)A*(p, -p - q)dp + c.c. 

h- AAn= f f ^MAO (p,q) f . 

8(27rf AypMl(p)Ml(q)M2(p + q)' ^"'^^ 

IssAAo{P:Ci) = |A(p, -q)|2 + A(p, -qfdpd^d^^^ + 2|A(p + q,p)p + 2A(p + q,p)'cipVp+q - 
4A(p, -q)A(p + q, p)dp - 4A(p, -q)A*(p + q, p)dqrf;+q + 
4A(p + q, p)A*(p + q, q)ciprfq + 4A(p + q, p)A(p + q, q)cip+q + c.c. 

c- AOO= ! ! IssAoojq.v) .25) 

8(27r)^yqypM20(p)M20(q)M2A(p + q)' ^ > 

IssAooiP, q) = |A(p, -q)P - A(p, -q)^cipVp+q + 2|A(p + q, p)P - 2A(p + q, p)=^(ipVp+q + 
4A(p, -q)A(p + q, p)dp - 4A(p, -q)A*(p + q, p)dqd;+q + 
4A(p + q, p)A*(p + q, q)dpdq - 4A(p + q, p)A(p + q, q)(ip+q + c.c. 

^ . nnn - / / Lsoooip. q) , . 

8 {2nf iq ip M20(p)M20(q)M20(p + q) ' ^'^^ 
Lsoooip, q) = |A(p, -q)|^ + A(p, -qfdpd^d^^^ + 

2A(p, -q)A(p, -p - q)dqO?p+q + 2A(p, -q)A*(p, -p - q)dp + c.c. 

2. Double bubble diagrams. The bubble diagram are shown in Fig.2e-g. 

e:AA = ^[ I ''^f'^^} (27) 
8(27r)VpyqMl(p)Ml(q)' ^ ' 

hbAAip, q) = A(-p + q, p + q)(ipfiq + 2/?p_q + c.c. 



/■ . -40 = ^ /" / -^bMo(P'q) ■ ^28) 
^ 8(27r)2iqipMM(p)APO(q)' ^^^^ 

46Ao(p, q) = -A(-p + q, P + q)cipciq + 2/?p_q + c.c. 
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8{27rf J^JpM^O{p)M^Oiq) 
hbOoiP, q) = A(-p + q, P + q)dpdq + 2/?p_q + c.c. 



3. Shift correction terms. The correction terms are given by 

8{2nfJJ,Ml{p)Ml{qy ^'"^ 
/cMA(p,q) = (2/?p- |7p|)(2/3q- |7q|). 

AO = ~^ f f ^^^-^o(P'q) ■ (o.) 

8(2vrf yjpMl(p)MS(q)' ^ 
/cMo(p,q) = (2/3p - |7p|)(2/?q+ |7q|)- 

8(27r)VjpMS(p)MS(q)' ^''^ 
^rOo(p, q) = (2/9p + |7pl)(2/9q + |7ql)- 



C. Cancellations of IR divergence within diagrams 

To analyze the IR divergence, one need to expand the propagators and vertices around 
small quasimomentum. Using the explicit expansion for A(ki,k2) given in Appendix A, 
one can in turn find those for 7k, (3^, and d^. It turns out that the constant and k'^ terms 
in M|(k) vanish, so that the (only) leading quartic term is — |a/i^|k|^ and MQ(k) = 



Oft 



2-|kp + (l-it)|kr 



. Here, 



P. 



St 



2 \a 



s+t 



^(^\ 5](-)^-rm*exp 



l.m 



{2^ 
2a2 



''A 



(33) 



As a result, the leading divergence ~ Jp Ipl^kt^q+pl"^ determined by the asymptotics 



of IssAAA{p,q} as both p and q approach zero. If hsAAA ~ 1, it would diverge as L^. 
However the vertex is "supersoft" at small quasimomenta so that the divergence is milder 
than expected. For example, the Aki^k;2^-ki-k2 vertex expansion at small momenta starts 
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with (/cf + /c^/cf ) . Therefore the leading divergence is "just" . Expanding IssAAa{P: q) 
around p = q = 0, we see that actually IssAAa{P: q) 0{p^) : 

IssAAAip, q) = Q/?00 - (322^ [pV " (p ' (p' " (p' + + 4p • 

As a matter of fact, it even becomes 0{p^°) after we symmetrize it with respect to p <-> q, 
and the diagram is actually finite. Similarly an apparent logarithmic divergence in setting 
sun AOO is nonexistent. 

D. Cancellation between different diagrams 

After the apriori most divergent diagram turned out to be convergent we look for milder 
IR divergencies in other diagrams. The remaining most divergent terms appear in contribu- 
tions ssaao, bbAA and ctaa, and come from the quasimomentum independent terms in the 
numerator of the integrands: 

hsAAO — 4/^00 ) hbAA — ^PoO, hrAA — Poo, 

respectively. Although they are divergent by themselves, their sum with appropriate 
coefficients — (2/5oo)~^, 1 and —1 cancels. The order divergencies come from the following 
integrands: 

IssAAO = -4/3o'o(p' + q^ + 2p-q), 

hbAA = Poo{-p^ -q^+p-q)- 

Expanding MQ(p-l-q) in the denominator to the second order in quasimomenta, we see they 
cancel each other after symmetrization with respect to p q and p — p. Finally, there 
are five are log(L) divergent terms: 

IssAAO = \Poo [9P' + + (P • g)(36p2 + 28q^) + 12pY + 36 (p • qf] - 
2/5oo/322 [2(p . q)(p' + q') - p'q' + 6 (p • qf] , 
hbAA = lf3oo [p' + q'- (2p • q){p' + q') + 6(p • qf] - 
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\l322 + - 4 (p ■ g) + g2) _ 6p2^2 ^ 2O ■ g)^] , 
hhAO = /^OO) 

hrAO = 3/5oo. 

By symmetrizing the sum of all the five terms, we see the final result is indeed free of IR 
divergence. 



IV. COMPARISON OF RESULTS WITH OTHER THEORIES AND 

EXPERIMENTS 

A. Comparison with high temperature expansion 

The same theory has been studied by various analytical and numerical methods. To 
compare our results with those obtained using other methods, let us restore the original 
units. The Gibbs free energy to two loops (finite parts of the integrals were calculated 
numerically) is 

^ = ~r — ^^'^ 
LzHhc 

9 = + ^ log(laTl) + c-^, (34) 

where numerical values of the coefficient is c = —5.2. Dimensionless entropy (the LLL scaled 
magnetization) is: 

dg f L.nc'mlb'Y^'' M 1 11.1 

- ' - ar - h ^c-g-, (35) 



and specific heat normalized to the mean field value is 

^ ^ 1 1 1 « 1 

+ 7^7r-6c-r- (36) 



j3A AC da\ (3a 27r 
We first compare the results with those of the high temperature expansion [H] . These series 



are known now to the 12th order in x Successive partial sums for specific heat at low 
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temperature are plotted on Fig. 3 (dashed lines) together with several orders of the high 
temperature expansion. Low temperature expansion indicates that the specific heat ratio 
grows with ax- On the other hand, the high temperature expansion clearly shows that 
it drops out fast above = 0. This means that there is a maximum in between which 
is consistent with most experiments and Monte Carlo simulations, see Fig. 4. Whether 
there is a melting phase transition either first order or continuous (in 2D it is necessarily 
of the Kosterlitz - Thouless variety) cannot be determined from series alone. Both series 
expansions have a finite radius of convergence, but this fact alone is not enough to decide 
that the singularity is at real temperature (it can be located in the complex plane as for 
example in the ID Ising model). The low temperature series are too short to estimate the 



radius of convergence. Naively comparing the second coefficient in specific heat eq.(|36D to 



the third one obtains an estimate = —yVhwl = —8. Phenomenologically first order 
melting occurs around = —10. Extensive analysis of the high temperature series has 
been made in |I^ , p^ ]. 



B. Comparison with MC, experiments and other theoretical results 

Low temperature results for free energy and magnetization agree well with available nu- 
merical simulations and experiments. However specific heat comparison is the most sensitive 
(second derivative). We therefore present here only the specific heat comparison. For 
< —5 the specific heat results on Fig. 4 are in accord qualitatively with experiments of 



T7| , and MC simulations of [jTO[- The same data were fitted by the theory of Tesanovic and 



coworkers IT^], We can calculate the coefficients of the low energy expansions from their 



theory and we get: 

-^ = --2- + 12/?A 
(3a AC Pa cit 



compared to our eq. (|36D . Even the sign of the l/a^ contribution is different. This theory is 
only an approximate one and perhaps some modifications are required in the low temperature 
limit. 
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V. UNDERSTANDING CANCELLATIONS OF INFRARED DIVERGENCIES IN 
A SIMPLE MODEL. NONPERTURBATIVE EFFECTS 



A. Toy model 

The dramatic cancellation of all the severe IR divergencies in the GL LLL model up 
to at the two loop level looks a bit mysterious. Although in critical phenomena the 
phenomenon of cancellation of " spurious" divergencies due to Goldstone bosons is well known 
p!8| , here it occurs under rather extreme circumstances. The theory is below its lower critical 



dimensionality. To better understand what is involved in these cancellations we investigated 
a model which has similar symmetry properties, but is much simpler. As was mentioned in 
the introduction the physics of the D dimensional GL theory is very reminiscent of that of 
the D — 1 dimensional scalar theory with two fields possessing a continuous 0(2) symmetry, 
see Fig. 5, 

where the dot denotes the derivative in the only dimension considered as "time" and we are 
interested in the spontaneously broken symmetry case aj- < 0. 

This model is equivalent to quantum mechanics of the two dimensional anharmonic 
oscillator. Of course one can solve this model nonperturbatively (albeit using numerical 
solution of the differential equation, we are not aware of the closed form of the ground state 
energy). Obviously the result is IR finite (bounded from below by the classical energy and 
from above by variational Gaussian energy) |ll9| . 



B. Perturbation theory 

It is important to trace the origin of the IR divergencies in the intermediate steps of the 
perturbative calculation. Although the calculation has been done using Feynman diagrams 
(steepest descent approximation of the path integral), like in GL model above, it is useful to 
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start from the standard time independent perturbation theory. Here we first have to decide 
what is the main part K and what will be a perturbation V . The perturbative vacuum in 
which < (pa >= Va is degenerate and we have many choices of the "unperturbed part. One 
of them corresponding to the choice 

Va = (V-or,0) = {v,0), 
O = ifi-v, A = (p2 

is: 

H^-f + K + V, 

V^2v{0^ + 0A^) + l{0^ + A'y 

where ttq and tta are conjugate momenta of the fields O and A respectively. The constant 
in H is the classical energy. Since any of these ground states are nonnormalizable, see Fig. 5, 
the IR cutoff L was introduced into K. It also removes the vacuum degeneracy. With this 
cutoff the unperturbed wave function is: 



'^o,om,<P2) = I ^ 1 exp 
Zero point energy (^o,o|^|^o,o) corresponds in the time dependent perturbation theory 



to 



Fi = ^ (TrlogGo + Tr logG^) ^ v + O (^^ 
in the time dependent formalism. The leading correction to the ground state energy is: 

It is equal to three "double bubble diagrams of Fig.2. The second order in V correction is: 

r,ss |(^0,0|^l^nl.n2)P 

-^2 - p ^ 

(nl,n2)^(0,0) 
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(^0,o|^|^l,o) P , I (^0,0|V^|^3,0) P , I (^0,o|^|^l,2) 



—2v —Qv —2v — I 

_ _3L2 LI _ 19 1 /I 

This contribution is the sum of the " setting sun diagrams" and correction terms (some terms 
which contain higher orders in l/ar were dropped). Unhke the GL theory there is no AAA 
setting sun diagram and therefore no L'^ divergence is expected. The leading divergence 
and the subleading L cancel between Fl'' and Fg***. 

C. Absence of long range order 

As is well known even discrete symmetry cannot be broken spontaneously in ID. This 
means that when we calculate perturbatively VEV of a quantity which is not invariant 
under the symmetry group 0(2) it should be IR divergent [Q. As an example we calculate 
expectation value of (pi. To first order correction to (\E'|(y9i|\E') arises from the corrected 
ground state 

^ = ^0,0 + ^ C„i,„2^nl,n2, 
(nl,n2)^(0,0) 

(^0,0|^|^nl,n2) 



Cnl,n2 



Cl,0 



2vL + 3 



5 



,3/2 



The result is 



1 L 



\/—aT 



It diverges linearly indicating "dynamical restoration" of the symmetry. The sign of the 
correction indicates that the VEV is reduced. The exact finite size expression is nonanalytic 
(hke a/— orexp ^J^z^ ) and approaches zero. 

This model clearly teaches us that although the 0(2) symmetry is unbroken, the per- 
turbation theory starting from the "broken" symmetry ground state not only cures its own 
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IR divergencies problems, but provides an accurate approximation to any 0(2) symmetric 
quantity. Perturbation theory actually "knows" about restoration of the symmetry missing 
only corrections of the essential singularity variety. 



VI. CONCLUSION 



To summarize, it is established up to order of two loops that perturbation theory around 
Abrikosov lattice is consistent. All the IR divergencies cancel due to soft interactions of the 
soft mode. Perturbative results as well as interpolation with the high temperature expansion 
agree very well with the direct MC simulation and experiments. The theory of Tesanovic et 
al [jlll has different low temperature limit and perhaps should be modified in this region. 

Consistency of the perturbation theory rules out a possibility of infinite lattice at any 
T > 0. Let us briefly summarize what it physically means and under what conditions this 
result is valid. The fact that v is IR divergent in both 2D (power) and 3D (logarithm) 
means that order parameter of translational symmetry breaking vanishes. Its correlator 
at very large distances approaches zero. This does not necessarily means that the state is 
liquid, namely the correlator decays exponentially with certain correlation length. It might 
decay only as a power like in 2D XY model ||2^ and melt into liquid via either first order or 
continuous transition. Assumptions are: no disorder, infinite sample, LLL approximations, 
no fluctuations of magnetic field. The third assumption can be relaxed since one obtains 



a supersoft spectrum also after including higher Landau levels, see ||l2|,|l3l, however 
including fluctuations of the magnetic field will probably stabilize the lattice in 3D since the 
spectrum becomes the usual Goldstone boson 1/k'^. 

Author is very grateful to B. Ya. Shapiro, T.K. Lee and D.P. Li for numerous discussions. 
Work was supported by grant NSC of Taiwan. 
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A. Appendix A. Small momentum expansion of Feynman diagrams 



In this Appendix we give formulas for expansion of integrands in powers of quasimomenta, 
which is needed to eighth order. The basic quantity which enter the Feynman diagrams is the 
function A. Using identities between the coefficients (3 defined in eq. (|33[) : /3o2 = |/5oo,/3o4 = 



j/3oo - 3/522, we have 



A(ki,k2)=exp^-Hli±M!^x 



Poo + \Poo [(^2 - ^kfY + {kl - iklf] + l/5oo [(^2^ - ik\)' + {kl - ikl)'] 

\P22 [-{kl ~ lk\f + {kl - iklf]' + 
o 

^(3o,[{kl-tklf + {kl-tklf] + 

^/524 {kl - iklf {kl - iklf [{kl - iklf + {kl - iklf] + 



(3,,[{k^^-^kyy + {kt-^kly 



40320 

^ P20 {kl - rklf {kl - iklf m - ^kl)' + {k't - iklf] + 



1440 
1 

576 
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FIGURES 
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FIG. 1. All the relevant propagators and vertices in the theory. Here, dash and solid lines 
represent the A and O fields, respectively. 
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FIG. 2. Setting sun and double bubble diagrams that contribute to the free energy at the two 
loop level. 




FIG. 3. Comparison between low temperature and high temperature expansion of the scaled 
specific heat defined in eq. 
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FIG. 4. Comparison between theoretical, experimental and Monte Carlo results of the scaled 
specific heat. 
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